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We reconsider a method based on dispersion theory, that allows one to extract the scattering 
length of any two-baryon system from corresponding final-state interactions in production reactions. 
The application of the method to baryon-baryon systems with strangeness S = —2 and S = —3 
systems is discussed. Theoretical uncertainties due to the presence of inelastic channels with near-by 
thresholds are examined for the specific situation of the reaction K~ d Tf^AA and the coupling of 
AA to the SN channel. The possibility to disentangle spin-triplet and spin-singlet scattering lengths 
by means of various polarization measurements is demonstrated for several production reactions in 
K~d and ^d scattering. Employing the method to available data on the AA invariant mass from 
the reaction ^^C{K~ , K'^AAX), a ^So scattering length of oaa = —1.2 ± 0.6 fm is deduced. 
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I. INTRODUCTION 



The baryon-baryon interaction in the strangeness S = —2 sector and, specifically, the A A system has been a topic 
of interest for quite some time. The fascination was generated not least by the possible existence of the so-called 
Q i/-dibaryon, a deeply bound 6-quark state with J — 0, isospin / = 0, and S = —2, predicted by R. Jaffe in 1977 
^ ' based on a bag- model calculation The binding energies of AA in nuclei, deduced from sparse information on 
doubly strange hypernuclei [^'^ indicated a strongly attractive ^5*0 A A interaction and seemed to be at least not 
inconsistent with the existence of such a bound state. The perspective changed drastically when in 2001 a new (and 
unambiguous) candidate for ^a-'^^ with a much lower binding energy was identified Q, the so-called Nagara event, 
^ ' suggesting that the AA interaction should be only moderately attractive. This conjecture concurs also with evidence 
, provided by various searches for the H dibaryon that did not yield any support for its existence, cf. Refs. for 
the latest experiments. 

However, very recently the iJ-dibaryon was put back on the agenda. Lattice QCD calculations by the NPLQCD 
[^, [l^ as well as by the HAL QCD [Tl|, [l^ collaborations provided evidence for a bound _ff-dibaryon. While the 
actual computations were performed for pion masses still significantly larger than the physical one, extrapolations 
suggest that the if-dibaryon could be still bound by around 0-7 MeV JJl] at the physical point, but it could also 
move above the AA threshold and dissolve into the continuum [l3, . (The original iJ-dibaryon [l[ was expected to 
be bound by roughly 80 MeV!) 

Indeed, the strength of the A A interaction as well as those of other S = —2 baryon-baryon systems is of rather 
k>( I general interest, noteably for a better understanding of the role played by the SU(3) flavor symmetry. Theoretical 
5_] ' investigations of the S = —2 sector that have been performed within the conventional meson-exchange picture [16142 ill , 
utilizing the constituent quark model ^2^ 2M : ^-l^o in the framework of chiral effective field theory (xEFT) [2J| 
all rely strongly on the SU(3) symmetry as guideline. Furthermore, the hyperon-hyperon (YY) interaction plays an 
important role in the understanding of the global properties of compact stars like neutron stars. Their stability and 
size as well as the cooling process depend sensitively on the strength of the YY interaction [2^ [2^ . 

As indicated above, practically the only experimental constraint we have so far on the AA interaction comes from 
the analyis of double-A hypernuclei. In the present paper we want to call attention to the fact that there is also another 
and even more direct way to determine the strength of the AA force but also the one in other S = —2 systems. It 
consists in studying the final-state interaction (FSI) of reactions where corresponding pairs of hyperons are produced. 
In fact, recently we proposed a method for extracting hadronic scattering lengths from production reactions (28l - 
[30| . The presentation of the method in those publications was done with special emphasis on its application to the 
hyperon-nucleon (YN) interaction. In particular, the reactions NN KYN and — KYN were analyzed, and 
possible uncertainties of the method were established. Polarization observables needed to disentangle different spin 
states of the final YN system were identified. 

In the present paper we explore the possibility of applying the method proposed in [28l - l30j to the AA system, 
but also other baryon-baryon states with S — —2 or even S = —3 are considered. Our study is motivated by the 
available data on the AA invariant mass distribution determined in the reaction ^^C(ii'^, X+AAA) [1, [3l|. These 
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data are afflicted by sizeable uncertainties, but still they allow us to demonstrate the practicability of our method and 
to extract an actual value for the AA ^Sq scattering length. In order to stimulate future dedicated experiments we 
consider specifically reactions like K~d — >■ KAA or K~d — >■ KEN where corresponding high-statistics measurements 
could be performed at J-PARC. The CLAS collaboration at Jlab has measured 7p — ^ K^K^E' fs?] and, thus, it 
might be feasible that they can perform also experiments for 'yd K^K^EN and 'yd -> K^K'-'AA Yet another 
option are reactions like pp — > K^K^AA and pp — > K'^K'-'Y.^J^^ which could be measured at the future FAIR facility, 
for example. Since the spin structure of these reactions differs partly from the ones considered in [28l-l30ll the question 
of what polarization observables are needed to disentangle the singlet- and triplet baryon-baryon states has to be 
re-addressed. Note, however, that for the considered AA system this issue is not relevant. Near threshold it can only 
be in the spin-singlet {^Sq) state. Due to the Pauli principle the other S'-wave, the ^Si, is forbidden in this case. 
Thus, no polarization experiment is required and, consequently, our method could be even applied to data on AA 
production on somewhat heavier nuclei, e.g. in K~ ^He or K~ ^He. 

Independently of that, the error estimation [2^ for the method has to be re-done. Specifically, for the AA system 
the inelastic threshold (due to the EN channel) lies with around 25 MeV much lower than for AA'' (where it is around 
80 MeV and due to the SiV channel). However, as we will see, the latter aspect increases the theoretical error of our 
method only marginally, under the discussed reasonable assumptions, and in the absence of a bound state. Taking 
into account the possibility of a bound state requires much more effort and is technically more complicated. In view 
of the current extrapolations of the lattice QCD results which rather seem to disfavor the existence of a bound state 
[l3l - [T5j we avoid the pertinent complications in the present study. 

For completeness let us mention that FSI effects as a tool to constrain the AA interaction were considered already 
many years ago by Afnan and collaborators [s^, H^- Their study was done under rather different presuppositions, 
namely for the reaction E^d ^ nAA and within the framework of Faddeev equations. With regard to the EN system 
there is also an entirely different possibility to determine the corresponding scattering lengths, namely via the study 
of E!~ atoms. Shifts of the energy levels due to the presence of the strong interaction would permit to deduce the 
scattering lengths for E~p or E~d, say, via the Deser-Trueman formula. The prospects of corresponding experiments 
were discussed in Ref. '3u\. 

The paper is structured as follows. In Sect. II we review briefly our method. Section III is devoted to the AA 
system. First we provide a new estimation for the error of the scattering length due to the extraction method, taking 
into account the relatively small separation of the AA and EN thresholds. Then we apply our method to available data 
on the AA invariant mass spectrum from a measurement of the reaction ^'^C{K~ , K^AAX). In Sect. IV we discuss 
several aspects of applying our method also to the EN and EE final-state interactions and even to the strangeness 
S = —3 sector. The paper ends with a short summary. Details of the polarization observables required for separating 
the spin singlet and triplet states are summarized in an appendix. 



II. REVIEW OF THE METHOD 



The basic idea of the method is to exploit the scale separation between a short-ranged production operator and a 
long-ranged final-state interaction (FSI). In this case the production operator can be regarded as point-like, and the 
FSI can be factored out. These conditions restrict the class of reactions and kinematic regimes that one can consider. 
Namely, one can only apply the method to reactions with large momentum transfer qt- Furthermore, the scattering 
length a in the system under consideration must have an appropriate magnitude, i.e. fulfil the condition a ^ l/<Zt- 
Sufficiently large scattering lengths are expected in the baryon-baryon sector. In particular, it is interesting to study 
the hyperon-nucleon and hyperon-hyperon interactions with different strangeness content of the hyperons. An elegant 
way to utilize the condition of scale separation is a dispersion-relation approach. Imposing unitarity and analyticity 
constraints on the amplitude and assuming that there are no bound states, one arrives at the following expression for 
the reaction amplitude As [281 1371,138} 
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where m is the invariant mass of the produced baryon-baryon system with the threshold value mo, s is the total 
center-of-mass (CM) energy squared, and t represents all the remaining kinematic variables the amplitude depends 
upon. The function ^(s, t, m'^) slowly varies with m^, which is a consequence of the assumed large momentum transfer. 
The cut off rrimax has to be determined in such a way that the integral extends over the whole region where FSI 
effects are expected to be important. Based on scale arguments a condition for nimax was derived in Ref. [2^ which 
reads, re-formulated in terms of the maximum kinetic energy in the two-baryon system, emax = 'nimax — 'm-o 
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Here as is the scattering length in question and fj, is the reduced mass of the baryon-baryon system. As argued 
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in Ref. '2q], for the hyperon-nucleon interaction a typical cut off is given by the condition Cmax ~ 40 MeV. The 
baryon-baryon scattering process has to be elastic in this region (i.e. there should be no other open channels) and it 
should be dominated by the s-wave amplitude parametrized by the phase shift Ssim?)- Note that formula ([T]) can 
only be applied to amplitudes for a specific baryon-baryon spin state S. Therefore, one has to be able to separate 
spin-singlet and spin-triplet states experimentally. The index "S" on the quantities above (and below) is a reminder 
that one has to consider the production of the baryon-baryon system in a definite spin S. 

It was shown in 28] how one can invert Eq. ^ to express the scattering length via the reaction amplitude squared 
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(or the differential cross section ^rf^^ ) 
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as = lim ^(^I^^)P r^'dm''- 

- ml (m'2 - m?) {p' (dm'^dt) }' 

where ma and mb are the masses of the two baryons, toq — ma + mi,, p' is the CM momentum in the baryon-baryon 
system and P indicates that the principal value of the intregral has to be taken. An analogous equation can be derived 
for the effective range. 

Possible theoretical uncertainties of the method originate from the following sources: (i) energy dependence of the 
production operator, (ii) influence of scattering at higher energies (m > mmax), (hi) contributions from inelastic 
channels (e.g. from the S7V o AN transition) and (iv) final state interaction among other pairs of particles. For the 
hyperon-nucleon FSI the theoretical uncertainty in the determination of the scattering length was estimated to be 0.3 
fm at most [5^. This estimate was confirmed by model calculations of production amplitudes using several different 
models for the hyperon-nucleon interactions with triplet and singlet scattering lengths varying from —0.7 to —2.5 fm. 

The general form of Eq. ([T|) admits approximations under certain conditions. One of the standard approximative 
treatments follows from the assumption that the phase shifts are given by the first two terms in the effective range 
expansion, 

pcomm^))^-- + ^-^p\ (3) 
a 2 

over the whole energy range, which is usually called the effective range approximation (ERA) . In this case the relevant 
integrals ([T]) can be evaluated in closed form as [s^ 

-l/a + (re/2)p2 - ip 



where a = l/?'e(l + — 2re/a). Because of its simplicity Eq. (|4|) is often used for the treatment of the FSI. A 
further simplification can be made if one assumes that a ^ Tg , a situation that is realized in the ^ 5*0 partial wave of 
the NN system. Then the energy dependence of the quantity in Eq. (|4]) is given by the energy dependence of the 
elastic amplitude 

A{m^) cx — — , (5) 

-l/a + {re/2)p^ - ip 

as long as p ^ l/?'e- Therefore one expects that, at least for small kinetic energies, NN elastic scattering and 
particle production reactions with a NN final state exhibit the same energy dependence [39l - l42| . which indeed was 
experimentally confirmed for meson production 'is'] . The treatment of FSI effects based on Eq. ^ is often referred 
to as Migdal- Watson approach [4^ |4l| , the one utilizing Eq. ^ as Jost-function approach. The reliability of such 
approximations as compared to the formula 1^ was investigated in detail in [2^. In general the method based on 
Eq. ([2]) works systematically better than the approximations and gives scattering lengths within 0.3 fm accuracy even 
for rather large scattering lengths like those for NN scattering. The uncertainty in the extraction employing the 
other two methods is typically larger. As demonstrated in [2^ these procedure lead to a systematic deviation from 
the true values of the scattering lengths of the order of 0.3 fm (Jost) and of 0.7 fm (Migdal- Watson) . 

III. THE AA SCATTERING LENGTH 

The AA system is certainly the most promising case where one could apply our method. Its threshold is the lowest 
one among all S = —2 channels and measurements could be performed for the reaction K~d — K'^AA, for example. 
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FIG. 1: Graphic representation of the thresholds for the various strangeness S — —2 channels with charges from Q — —2 to 
Q = +2. The thresholds for isospin averaged masses are indicated at the bottom. 



Moreover, no polarization experiment is required because (near threshold) the AA can be only in the (spin singlet) 
^Sq partial wave. The '^S'l state is forbidden due to the Pauli principle, as already mentioned. Thus, spin triplet 
states can only occur in P (or higher partial) waves - and it is save to assume that such higher partial waves do not 
contribute near threshold. There is, however, a complication because the first inelastic threshold (due to EN) is fairly 
close: the E^n channel opens at an excess energy of 23.06 MeV, c.f. Fig.[TJ (In the Ap case considered in [28| the first 
inelastic channel (SV) opens at 76.96 MeV!) Thus, it is necessary to re-do the error estimate of Ref. [2^. This will 
be done in subsection A below. Note that, for convenience, we will work with isospin-averaged masses throughout 
this section so that the SN threshold is located at 25.8 MeV! 

In subsection B we apply our method to the AA invariant mass distribution measured in the reaction 
'^'^C{K-,K+ AAX) by the KEK-PS E224 Collaboration [8|]. Those data, though afflicted by sizeable error bars, 
allow us to demonstrate how our method works, and they even enable us to deduce a concrete value for the AA ^So 
scattering length. 



A. Error estimation 



In this subsection, we generalize the discussion of theoretical errors presented in [23| to the case of the occurrence of 
inelastic channels. Specifically, we estimate the theoretical error for the extraction of the hyperon-hyperon scattering 
length exemplary for the reaction K~d — >■ K^AA taking into account that there is a near-by inelastic threshold due 
to the coupling of AA to the 57V channel. 



The uncertainties originate 28| from the energy (i.e. m^) dependence of the function ^{s,t,'m^) in Eq. ([T]). These 
include the energy dependence of the production operator (i.e. the influence of left-hand singularities), contributions 
of the elastic scattering to the dispersion integral at higher energies, the influence of inelastic channels, and the 
interaction between other pairs of particles in the final state. The latter effect can be controlled by choosing different 
kinematical conditions such as initial energy (flnal-state interaction among other pairs of particles would depend on 
such a choice whereas the A A FSI does not). Also investigating the invariant- mass distribution for a corresponding 
pair of particles via a Dalitz plot analysis can provide additional information on their interaction p3l |. In what follows, 
we will disregard this (possibly important) kind of correction and focus on the other three. 

The energy dependence of the production amplitude A5(s,t,m^) can be deduced from the basic principles of 
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analyticity and unitarity. The discontinuity of the amplitude in m is given by the sum of the elastic term, the 

?™(s, t, m?), and the left-hand part, D^g^ 



inelastic contribution of the SA^ (and/or other) channel, D™{s, t, m^), and the left-hand part, Dg'^-{s, t, m^), denoting 



the remaining contribution from the production operator 
Ds{s,t,m^) = -^{Asis,t,m^ +iO) - Asis,t,m^ - iO)) = A{s, t, m^)e-'^ sinS + D^g''-{s,t,m'^) + Dp{s,t,m^). (6) 
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The inelastic contribution reads 



/ A(s,t,m^)(l-77)e-^^'^ 



2i 



P2 



m"^ — (in non-relativistic kinematics), 



(7) 



where r] is the inelasticity parameter in the AA system, A2 is the production amplitude for the reaction K~d — )■ K~^EN, 
and /12 is the AA — >■ transition amplitude. The latter can be written in terms of 77, (5, and the phase shift 62 of 
the EN channel: 



/l2 = 



2^/Rp2 



(8) 



In order to shorten the notation we rewrite as 



1/12^)1 ^-.j 



/i2(m2) = 2p2/i2(to'), ~S^S + S2+Sa-Sa, 



(9) 



where we suppressed any dependence on s and t, and where we denoted the phases of the production amplitudes A 
and Ao by ^4 and 6a2, respectively. The solution of Eq. ([6]) in the physical region can be represented as (see Refs. 
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where rfi^ denotes the upper end of the left-hand cut. In order to remove the energy independent part from the inelastic 
dispersion integral, we made subtractions at to^ — itiq in the definition of u(m^) and in the inelastic dispersion integral 
(a constant term is assigned then to the left-hand contribution which anyway, as we will show, is slowly varying with 
energy.) 

The theoretical error of the extracted scattering length is determined by the energy dependence of the function 
$(s,i,TO2) from Eq. ^ [H] 



lim P 
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(11) 



The function $(m^) depends on m^^^^. This dependence factors out $(w^j,^ 
\I'(m^^^,m^) contains the information on the phase shift at energies above mmax 

<5(to'2) 



^fm^ ™2 



TO ) oc exp 



2 — TO^ — iQ 
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^'(TOj^jj^, to^)$(to^) where 



(12) 



Thus the theoretical error is the sum Sa'^^^^ — Ja™™°^ -|- 5d, where (5a™'""=' is due to the factor ^'(to^q^, to^) and 5a is 
determined by the energy dependence of $(to^) which is related to the energy dependence of the production operator 
and inelastic effects. 

Let us first estimate the theoretical error originating from the energy dependence of the production operator, 
neglecting for the moment the inelastic contributions. In order to do this one has to investigate the contribution to 
the dispersion integral from the left-hand cuts. We follow here the procedure utilized in [28|. We are interested in 
the left-hand singularities of the amplitude, i.e. singularities in some momentum transfer variable t. The simplest 
production mechanism for the reaction K~d — >■ K'^KK is seemingly the one shown in Fig. [5] denoting the exchange 
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FIG. 2: Typical production mechanism for the reaction K d ^ K AA. 



of one nucleon and one pion (one should add, of course, the diagram with proton and nucleon interchanged and 
7r+ replaced by tt^). Clearly, there are also other more complicated production mechanisms that will contribute. 
However, those should be of even shorter range and thus, correspond to production operators with even weaker 
energy dependence so that they are not relevant for the estimation of the theoretical error. We consider for simplicity 
kinematics corresponding to the final particles produced not far away from threshold (so that the momenta of the 
final particles are small compared to the momentum transfer) . The required initial CM momentum is then pi ft! 580 
MeV/c. The production operator contains a cut corresponding to the tt~p intermediate state, which can be associated 
with the interaction in the AK'^ system discussed above. The other cut over the neutron and tt^ is the one we are 
interested in. We can very roughly estimate the energy dependence associated with this singularity via approximating 
it by a pole term, 

$(^2) _ t = {pk- -Pk? , (13) 

t — _ 

n-TT 

where pi^- and pA are the corresponding 4-momenta and m„^- is the effective invariant mass of the n-K~ system, 
i.e. rrin-K- ~ 1 GeV. For threshold kinematics (i.e. zero momentum of all final particles - this choice is made for 
definiteness) this energy dependence (after averaging over all directions) has the form $ ft 1 + p'^/p'^, where p is the 
CM momentum in the AA system and pt is of the order of 2 GeV. The correction to the scattering length due to 
such an energy dependence of the amplitude amounts to [1^ ^^f^ ft 0.01 fm for values of emax = 40 — 60 MeV. The 
above rough estimation is sufficient to observe that the error coming from the energy dependence of the production 
operator is negligible as compared to the other sources of uncertainties and we can savely assume that 

$/./..(m2) ft ^i.hXml) = A{ml) ft A(m2)e-"("^'). (14) 

Note that since the AA scattering length is expected to be somewhat smaller than OAp E^, H^, [13] we choose 
^max = 60 MeV rather than emax — 40 MeV (used for AA^ scattering) as our central value in order to minimize the 
effect of higher energy scattering. As was pointed out in Ref. [2^ Cmax must be chosen well above (cf. also 

the discussion in Sect. |ll|. Although a further increase of emax can help even more in reducing this effect, in reality 
it would be difficult in this case to separate S'-waves from higher partial waves in the final state and to avoid the 
influence of the interaction in other channels - the effects and the number of inelastic channels would increase. 

Next we consider the error coming from the inelastic channels coupled to the AA system. For simplicity we consider 
only the one nearest to the AA threshold, namely E,N, which opens only about 25 MeV above the AA threshold so that 
it is necessary to analyze its impact. Clearly, due to the lack of empirical information on the 'E.N interaction and the 
coupling of this channel to the AA system such an error analysis cannot be done in a completely mo del- independent 
way. One has to make some assumptions on the strength of the interactions in the relevant channels. For our analysis 
we prepared three variants of the hyperon-hyperon interaction of [2^ which yield AA scattering lengths of —1.36 fm, 
— 1.50 fm, and —1.70 fm, respectively. We are interested in the situation when the effect coming from the inelastic 
channel is small and its contribution can be treated perturbatively, i.e. 

(m 2) ("^')- (15) 

If this is not the case the error of the extraction of the scattering length would be comparable with the scattering 
length itself and then an extraction would no longer be meaningful. We made a subtraction in Eq. (|10p at rri^ = uiq so 
that is exactly zero at the beginning of the integration interval. Then a small resulting correction to the scattering 
length would imply a likewise small variation of so that then the condition (|15p would be justified. Note that 
formally the subtraction point does not enter the expression for the error. 

Utilizing Eqs. and one can rewrite the formula ([TT|) for the inelastic contribution to the theoretical error 
in the form 
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Substituting ^^(to^) from Eq. (jlOp and performing one integration one gets 
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where Prnax = 



Using the exphcit form of D'^gim?) from Eq. (O and making a suitable change of variables we have 
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From [2^ we recall the expression for Ja™""^ (with 
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Note that the integration over y from to cx) can be truncated at y = 1, say, for practical applications since the 
integrals are rapidly converging (unless the phase shift is rising unnaturally fast with energy) . In any case, we cannot 
trust the ^EFT predictions for the amplitude at such high energies. 



The numerical values of ,5 and of | ^ | are not known, therefore we estimate 
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where Max | ^ | is the maximal value of this ratio in the considered energy region. 

Using the three mentioned variants of the hyperon-hyperon interaction we arrive at the following estimates of the 
theoretical errors for emax = 60 MeV: 
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A. 



A. 



A. 



A 



f m . 



f m . 



f m . 



8 




FIG. 3; Dependence of the extracted scattering lengths on the value of the upper limit of integration, e,nax- Shown is the 
difference to the exact results for three variants of the xEFT AA interaction with a — —1.50 fm (solid line), a = —1.70 fm 
(dashed line), and a = —1.36 fm (dotted line). The shaded area indicates the estimated error of the applied method. 



for the variant with a = —1.70 fm. The value of Max |^| can only be accessed from the corresponding production 
experiment for the 'EN channel. Some estimates can be obtained by looking at the strength of the cusp effect in 
the AA production channel. Under the assumption that there is no specific production mechanism that makes this 
ratio large, one can estimate the ratio, at least qualitatively, from the unitarity contribution that correspond to the 
(on shell) AA — >■ SiV conversion. It amounts to Max|^| = Max|/i2P2| < 0.3 for all three considered interactions. 
Based on those numbers a rough estimation for the full theoretical error related to inelastic effects and higher energy 
scattering yields 6a*^ < 0.3 — 0.4 fm. Such small values for (5a*'* justify the approximation made in Eq. because 
it means that $i„(TO^) does not change much for p < an energy range that savely covers the region we are 

interested in. 

In order to test our method and to check our error analysis we applied it to a production amplitude, calculated 
with a point-like production operator, that incorporates AA final-state interactions generated from the three variants 
of the xEFT interaction. In Fig. [3] we show the dependence of the difference of the extracted scattering length and 
the exact one on the cut off emax of the integration. One can see that for Emax = 60 MeV the theoretical error is 
indeed within the range of 0.25 — 0.35 fm in agreement with the preceding analysis. Note the apparent drop of the 
curves around 25 MeV, i.e. at the opening of the EN channel. 



B. Analysis of data on the AA invariant mass from ^'^C{K-,K+AAX) 

First results for the AA invariant mass distribution were reported by the KEK-PS E224 Collaboration from a 
measurement of the ^'^C{K^ , K^AAX) reaction in 1998 |3l|. An enhancement was seen for invariant masses near 
threshold. Already at that time there were attempts to extract the AA interaction from the spectrum [4^, |4^. Then 
in 2007 the KEK-PS E522 Collaboration published a AA invariant mass spectrum with somewhat better statistics 
. Also in this case efforts were made to extract the AA scattering length. The value reported at some conferences 
4Y|, |48| , CAA = -0.10t?;^^±0.28 fm, was obtained by utilizing the Migdal- Watson approach (Eq. ([5])). 

As already pointed out above and as we thoroughly investigated in [29|], the Migdal- Watson approach works only 
well for fairly large scattering lengths, i.e. for values of the order of 5 fm or more, as they are typical for the NN 
interaction. For small scattering lengths as suggested by the analyis in [U 13 ^^^^ approach is not reliable. It can 
lead to a systematic deviation of 0.7 fm or more. Thus, we re-analyzed the AA invariant mass spectrum given in 
3, employing our method based on Eq. Indeed, since in the AA case the FSI can only occur in the ^5*0 partial 
wave and, thus, no polarization experiment is required, our method can be applied also to data like those of AA 
production on carbon Q from the reaction ^^C{K~ , K^AAX). But one has to keep in mind that in reactions on 
nuclei the energy dependence of the production operator is not so well under control. For example, there could be 
excitations in the other fragments of the reaction process. In addition, in the concrete case, the error bars of the AA 
invariant mass distribution are quite large, therefore one has to expect large uncertainties for the extracted scattering 
length. Nevertheless, in order to demonstrate how the method works we applied it to the data of Ref. Q, following 
the procedure for the analysis of experimental data described in detail in the Appendix A of We fit the data 
with the amplitude squared parametrized as 



|^(m)p = exp 



__ci 

(m2 - C'i 



Co + 7 9 



(21) 
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m^-2m^[MeV] 



FIG. 4: AA invariant mass spectrum for the reaction ^^C{K , K'^ AAX) 8] and our fit to it (solid line) 



multiplied with the phase-space factor, and allowing for a finite mass resolution of 2.5 MeV. The resulting curve is 
shown in Fig. |4] (solid line). Then we use this fit to extract the scattering length from the dispersion integral with the 
cut off emax — 60 McV. The result is 

AAA = -1.2 ±0.6 ±0.4 fm, (22) 

where the first error is due to the uncertainties in the data and the second value is the theoretical error estimated in 
the preceeding subsection. 



TABLE I: AA ^So scattering lengths (as) and effective range parameters (rs) for various strangeness S = —2 interaction 
potentials (in fm). In case of the xEFT interaction results for the lowest (550 MeV) and highest (700 MeV) cut-off value are 
given, cf. (2^ . Note that (a) the scattering lengths of the Nijmegen (ESC04) potential differ significantly depending on whether 
they are calculated in particle ^ITJ or isospin [18j basis, (b) in the potentials by Tominaga et al. ^201 some channel couplings 
are not included. 



YY interaction 


reference 


flj [fm] 


rs [fm] 


XEFT (550) 


[24] 


-1.52 


0.82 


xEFT (700) 




-1.67 


0.34 


Nijmegen (NSC97a) 


[16] 


-0.27 


15.00 


Nijmegen (NSC97f) 


[16] 


-0.35 


14.68 


Nijmegen (ESC04a) 


m 


-3.804 


2.42 


Nijmegen (ESC04d) 


[18] 


-1.555 


3.62 


Nijmegen (ESCOSa") 


[19] 


-0.88 


4.34 


Tominaga (set B) 


[20] 


-3.40 


2.79 


Fujiwara (fss2) 


m 


-0.821 


3.78 


Valcarce 


[23] 


-2.54 





For the ease of comparison we present in Table H] a selection of AA scattering lengths predicted by various 
published YY interaction models. The large scattering length of the interaction by Tominaga et al. from 1998 [l^ 
was still triggered by the first experimental information on the ground states of AA^e, aa^c, and aa^ [1-13 ■ Those 
experiments suggested a separation energy of ABaa = 4 — 5 MeV, where the separation energy is defined as, e.g., 
Ai?AA = -BAA(AAHe) — 2i?A(AHe). Such a large separation energy could be only described with a rather strong AA 
interaction with scattering lengths caa in the order of -2 to -3.6 fm [2O, [4^ or larger [50| . 

The analysis of the unambiguously identified aa-Hg double hypernucleus (Nagara event) [Ej yielded the much smaller 
separation energy of 1.01 ± 0.20 MeV. Calculations that obtain separation energies in agreement with the new exper- 
imental value suggest AA scattering lengths in the order of —0.7 to —1.3 fm [itI. [2^ ISOj. It should be said, however, 
that so far there are no fully microscopic (i.e. six-body) calculations of aa-^g available that utilize only elementary 
baryon-baryon {NN, YN, YY) interactions such as those listed in Table ID All studies are performed either with 
three-body Faddeev equations applied to the cluster model [scMs^l , the Brueckner theory approach [13, [E^ , or with 
the stochastical variational method 56], and rely, at least partly, on effective two-body interactions. 
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In this context let us mention that the separation energy for the Nagara event has been recently re-analysed and is 
now given as ABaa = 0.67 ± 0.17 MeV [57]. See also Ref. f58| for a recent review of the status of AA hypernuclei. 

TABLE II: EE and HA'' S-wave scattering lengths a and effective range parameters r for various strangeness 5* = —2 interaction 
potentials (in fm). The subscripts s and t refer to the singlet ( ^So) and triplet ( ^Si) states, respectively. In case of the xEFT 
interaction results for the lowest (550 MeV) and highest (700 MeV) cut-off value are given, cf. [24 ]. Note that (a) the scattering 
lengths of the Nijmegen (ESC04) potential differ significantly depending on whether they are calculated in particle [itI ] or 
isospin [3 basis, (b) in the potentials [20l . [2H some channel couplings are not included. 



YY interaction 


reference 


channel 


as [fm] 


Ts [fm] 


at [fm] 


rt [fm] 


XEFT (550) 


[24] 


E+E+ 




-6.23 


2.17 


- 


- 


XEFT (700) 


m 


E+E+ 




-9.27 


1.88 


- 


- 


Nijmegen (NSC97a) 


[16] 


E+E+ 




10.32 


1.60 


- 


- 


Nijmegen (NSC97f) 




E+E+ 




6.98 


1.46 


- 


- 


Fujiwara (fss2) 


[22] 


E+E+ 




-85.3 


2.34 


- 


- 


Valcarce 


[23] 


E+E+ 




0.523 




- 


- 


XEFT (550) 


[24] 


E°n 




- 


- 


-0.34 


-5.86 


xEFT (700) 


]2£ 






- 


- 


-0.15 


16.3 


Nijmegen (ESC04a) 


\m 


EN{I = 


0) 


- 


- 


-1.672 


2.70 


Nijmegen (ESC04d) 


[18j 


HA (J = 


0) 


- 


- 


122.5 


2.083 


Nijmegen (ESCOSa") 


m 


HA(7 = 


0) 


- 


- 


6.9 


1.18 


Tominaga (set B) 


[20] 


HA (J = 


0) 


- 


- 


-0.352 


17.4 


Ehime (1.82) 


[21] 


HA(/ = 


0) 






-0.43 


13.0 


Valcarce 


[23] 


HA(/ = 


0) 






0.28 




XEFT (550) 


[24] 


H^n 




0.21 


-30.7 


0.02 


968 


XEFT (700) 


[24] 






0.13 


-98.5 


0.03 


548 


Nijmegen (NSC97a) 


[16] 


H~n 




0.46 


-6.09 


-0.04 


634 


Nijmegen (NSC97f) 


[16] 


H^n 




0.40 


-8.88 


-0.31 


870 


Nijmegen (ESC04a) 


[18] 


3~ n 




0.491 


-0.421 






Nijmegen (ESC04d) 


[18] 


H~n 




0.144 


4.670 






Nijmegen (ESCOSa") 


fl9] 


H^n 




0.58 


-2.71 


3.49 


0.60 


Tominaga (set B) 


[20] 


H~n 




-0.202 


33.0 


-0.484 


10.6 


Ehime (1.82) 


[21] 


H~n 




-0.27 


20.3 


-0.56 


9.0 


Fujiwara (fss2) 


[22] 


H^n 




0.324 


-8.93 


-0.207 


26.2 


Valcarce 


[23] 


H«p 




-3.32 




18.69 





IV. THE HA AND EE SCATTERING LENGTHS 

Since the AA system is a pure isospin / = state, the 1=1 EN interaction is also elastic and, thus, permits a 
determination of the corresponding scattering length via our method. The first inelastic channel, AE, opens at the 
excess energy of around 52 MeV and, therefore, should affect the extraction of the scattering length less than what 
has been discussed in the context of the AA case above. The required EN invariant mass spectrum is accessible 
experimentally in the reaction K~d — >■ K^E~n, for example. However, since EN can occur in the ^Sq as well as 
in the '^Si partial wave, one needs data from an experiment with polarization for a separation of the singlet and 
triplet contributions. The relevant observables are discussed in Appendix |Al In addition, one has to keep in mind 
that the reaction K^d K^EN might be dominated by the quasi-elastic process K~N — K~^E, similar to what 
happened for the reaction KYN [s^] with ^N — KY . In such a case the reaction kinematics has to be chosen 

rather carefully in order to suppress the contributions from the quasi-elastic process, as studied in detail in [30l [. which 
certainly increases the difficulties for a corresponding experiment. 

As a subtlety let us mention that even in the E^n system the scattering length for '^Si is real, although the amplitude 
is actually the sum of / = and 1 = 1 states. Because of the Pauli principle, the '^Si partial wave of the AA system 
is forbidden so that there is no coupling between the AA and EN channels for the partial wave in question. 
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Experimental information on the 2iV invariant mass spectrum is scarce |59l - l6l| . Results published in Ref. [59[ 
for the case, obtained from a K~d bubble-chamber experiment, suggest an enhancement in the invariant mass 
distribution - but at around 2480 MeV and not near the EN threshold which is at around 2255 MeV. In any case, the 
statistics is too low for drawing any conclusions. The situation is better for an experiment performed at Saclay jg^] 
where the missing mass (MM) in the reaction K^d — > + MM at 1.4 GeV/c was studied. The curve presented in 
this publication exhibits a rather smooth behaviour around the SA^ threshold which suggests that the S^n interaction 
in the ^Sq and/or ■^S'l might be fairly weak. Such a conjecture is actually in line with the results of several of the 
potential models summarized in Table HIl which predict rather small scattering lengths. We want to emphasize, 
however, that there are also models with a fairly strong S7V interaction. Noteably the latest Nijmegen potential 
(ESCOSa") produces bound states in the '^S'l partial wave of the / = and 1 = 1 channels [l^. The binding energies 
are comparable to that of the deuteron and, accordingly, sizeable near-threshold enhancements in the corresponding 
invariant mass spectrum are to be expected if such bound states indeed exist in nature. 

Besides the E.N system the 1 = 2 channel is potentially interesting too because it is also elastic. But due 

to the charge it cannot be produced with a beam on the deuteron. However, the S^S^ system which is likewise 
1 = 2 could be studied, namely in the reaction K~d — K^tt^Y,~Y,^ . Also in this case, only the ^Sq partial wave is 
present so that no polarization data are required for a determination of the scattering length. Contrary to S^n, here 
practically all model predictions for the (E^E^) scattering length are fairly large, cf. Table HIl Details of the 

application of our method to cases where the Coulomb interaction is present can be found in Ref. 29] . 

In principle, one can even consider reactions of the type K~d — >■ KK'BA and K^d — >■ KK'E.Yi, which would give 
access to the strangeness S = —3 world. Potential-model calculations [U [13 and also predictions obtained within 
the framework of xEFT [62] suggest that the interaction in some of the channels are strongly attractive so that the 
corresponding scattering lengths could be large. Clearly, here it would be desirable to have at least a rough estimation 
of the count rates that one can expect in order to judge the feasibility of such experiments. Independently of that, 
also for these reactions we provide and discuss the relevant polarization observables that are needed for a separation 
of the singlet and triplet contributions, cf. Appendix \K\ 



V. SUMMARY 



We reviewed a method that allows one to extract hadronic scattering lengths from production reactions by studying 
final-state interactions. In particular, we discussed its applicability to the case of baryon-baryon interactions in the 
strangeness S = —2 and S = —3 sectors. We emphasized the importance of separating different spin states of 
the interacting particles. Considering as examples the reactions K^d — > KB1B2, jd K1K2B1B2 and K~d — 
K1K2B1B2, we could demonstrate that it is possible to construct polarization observables that provide access to 
spin-singlet and spin-triplet scattering lengths. In case of the AA and (or S~S~) interactions near threshold 

only the ^Sq partial wave is present due to the Pauli principle and, thus, no polarization experiments are required for 
determining the pertinent scattering length from the final-state interaction. Employing the method to available data 
on the AA invariant mass from ^^C{K^ , K'^AAX) a ^5*0 scattering length of oaa = —1.2 ± 0.6 fm is deduced. 
The error given here reflects the accuracy of those data. Thus, it would be important to perform experiments with 
better statistics which could be done, e.g., at J- PARC ^635. This would then allow one to reduce the error on the AA 
scattering length to the one of the extraction method, which we estimate to be in the order of 0.3 — 0.4 fm. 
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Appendix A: Spin considerations for the production of baryon-baryon systems with strangeness S = 2 and 

5" = -3 

The technique utilized for the hyperon-nucleon interactions in [28l - [30| is applicable also in the case of strangeness 
S = —2 and S = —3 systems. The necessary condition that baryon-baryon scattering should be elastic up to some 
TO = rumax IS Satisfied for the following S = —2 channels: AA, E~E~, S^n, and for the S = —3 channels: 

S~A, S^A, S~E~, S'^E+. See Fig.[T]for a graphic overview of the kinematics. In what follows we are going to consider 
as examples K" d, pp, and jd scattering in complete analogy with the hyperon-nucleon production reactions studied 
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in [30j. The following four types of reactions can be used to produce baryon-baryon states with strangeness S = —2 
and S = —3: 

K-d^KBiB2 {K"AA, K+E-n), (Al) 
pp^KiK2BiB2 {K+K+ AA, K+K°E°p, K°K"J:+Y.+), (A2) 
-fd K1K2B1B2 {K+ K° AA, K°K°E°p, K+K+E-n), (A3) 
K-d^KiK2BiB2 {K°K°E"A, K+K"E-A, K+K+E-Y.-). (A4) 

Now we come to the question of separating the spin-triplet and spin-singlet states in the baryon-baryon system. As 
in the case of the hyperon-nucleon interaction it is sufficient to consider reactions with polarized initial particles. We 
start from the general form for the reaction amplitude in the center-of-mass (CM) system for the three processes 

MK-d^KB^B, = alicd X p) ■ k + a\{ed ■ S') + al{ed-p){S' ■ k) + a\{ed ■ k){S' ■ p) , 



Mpp^KiK2BiB2 = &i + blip X fc) • 5 -I- blip ■ k){p xk)- S' + {h\p,k.j + h\pjh)S^S'j 



Myd^KiK2BiB2 = cl{ey X Ed) ■ p + C2{e-y ■ ed){S' ■ p) + cl{e^ ■ S'){ed ■ p) , 
MK-d^K,K2B^B2 ^ d{{ed-p)+d\{edxS')-p. (A5) 

Here a, 6, c, and d are some functions of s and of m (and, in general, of further invariants that are required to 
specify the kinematics of the reaction), where their upper indices indicate whether they correspond to spin-singlet 
(s) or spin-triplet {t) amplitude. The polarization vectors of the deuteron and photon are denoted by td and e!y, 
respectively. The spin vectors 5, S' are used for the spin-triplet initial and final states, respectively. For the last 
two reactions we assume the momenta of the final kaons to be either aligned or anti-aligned with the direction of the 
initial CM momentum p. This leads to a significant simplification allowing one to separate different spin states. For 
the reaction K~ d — KB1B2 such a restriction is not necessary and the momentum of the emitted kaon is denoted 
by k. For the reaction pp — > K1K2B1B2 we assume for simplicity that both emitted kaons go into the same direction 
k. It is convenient to introduce the following set of polarization observables 

Oi = (1 - ^/2n)3^ , ^ (2 + V2T0„)3^ , O3 = T{1 



'dm^dt' ^ ' '"'dm^dt' ^ ^"dm^df 

O4 = h + V2T^o + V3in2+Tl2)):^=V3{-V2T^, + {n2 + Tl2)) 



where the various T's for the ^d initial state are defined in |30|, and -j^r^ is the unpolarized differential cross 
section. T^q and T^q have the same definition also for the K~d initial state, the only difference being the absence 
of the summation over the photon polarizations. The observable Oi selects the amplitudes with longitudinal target 
polarization ed \\ p, whereas O2, O3, O4 select the amplitudes with _L p. In addition O4 contain only that part 
of the amplitude which is antisymmetric with respect to an interchange of ed and e^. The observables O5 and Og 
correspond to the proton-proton induced reaction. Here Aoy is the analyzing power and Ayy is the spin correlation 
coefficient for polarized beam and target (2S . i43j , and y is the direction perpendicular to the reaction plane. 

Now inspecting the structure of the reaction amplitudes (jA5l) we can identify the observables that allows one to 
separate a particular spin state: For the reaction K~d — > KB1B2 the triplet final state can be singled out by the 
observable Oi for any direction of the emitted kaon, or one can measure the unpolarized differential cross section 
for k II p. The spin singlet state cannot be separated. For the reaction ^d — >■ K1K2B1B2 the triplet final state can 
be separated by measuring Oi. The observable O4 provides access to the spin-singlet amplitude. For the reaction 
K~d — )■ K1K2B1B2 the observable Oi separates spin-singlet contribution, whereas O2 and O3 separate spin-triplet 
state. For the reaction pp — K1K2B1B2 the observable Oq separates the spin-triplet contribution, whereas O5 is 
proportional to sin0 3{&f&2* + ^4^5* cosfl} with cos9 = p ■ k. Since the bi's are even functions of cos9 (due to parity 
conservation), after the integration of O5 over an angular region symmetric with respect to 6 — ^ only spin-singlet 
amplitudes survive. 
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